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Discovery of 172-s Pulsations from a Be/X-ray Binary Candidate 
AX J0051.6-7311 in the SMC with ASCA 



O 
O 

o 

(N 

Or 

<D 



> 

o 

On 
O 
O 
O 



Oh' 
i 

O 

CO 

c3 



X 



Jun YOKOGAWA, 1 Ken'ichi Torii, 2 Kensuke Imanishi, 1 and Katsuji Koyama 1 * 

1 Department of Physics, Graduate School of Science, Kyoto University, Sakyo-ku, Kyoto, 606-8502 

E-mail(JY): jun@cr.scphys.kyoto-u.ac.jp 
2 National Space Development Agency of Japan, 2-1-1 Sengen, Tsukuba, Ibaraki, 305-8505 

(Received 2000 July 17; accepted 2000 August 8) 
Abstract 

The results from three ASCA observations of AX J0051. 6-7311 = RX J0051.9-7311 are reported. 
Coherent pulsations with a barycentric period of 172.40 ± 0.03 s were discovered in the third observation, 
with an exceptionally long exposure time of ~ 177 ks. The X-ray spectrum was found to remain unchanged 
through these observations, with a photon index of ~ 0.9. Energy-resolved pulse profiles in the third 
observation reveal that the pulsations are mostly due to photons with an energy above ~ 2 keV. Archival 
data of ROSAT and Einstein indicate that AX J0051.6— 7311 exhibits a flux variation having a factor > 20. 

Key words: pulsars: individual (AX J0051.6— 7311) — stars: emission-line, Be — stars: neutron — 
X-rays: stars 



1. Introduction 

Candidates for X-ray binary pulsars (XBPs) in the 
Small Magellanic Cloud (SMC) have been identified in 
two ways. Yokogawa et al. (2000a) observed the main 
body and the eastern wing of the SMC with ASCA, 
and made systematic analyses on the detected X-ray 
sources. They found that XBPs and thermal supernova 
remnants are clearly separated by the spectral hardness 
ratio. Eight candidates for XBPs were found from this 
classification. Probably, no detection of coherent pulsa- 
tions from these candidates is due to the limited photon 
statistics. On the other hand, Cowley et al. (1997) and 
Schmidtke et al. (1999) made a systematic search of op- 
tical counterparts for X-ray sources. They determined 
accurate X-ray positions by ROSAT/HRI and made op- 
tical photometry and spectroscopy within the X-ray error 
circles. They identified five X-ray sources with a Be star 
companion, which are good candidates for XBPs. 

The most promising and direct way to establish an 
XBP is to search for coherent pulsations from these can- 
didates. ASCA observations are suitable for a pulsation 
search because of high sensitivity in the hard X-ray band. 
In fact, we searched for coherent pulsations with ASCA 
from RX J0050.8— 7316, a binary with a Be star com- 
panion (Cowley et al. 1997), and discovered 323-s pulsa- 
tions (Imanishi et al. 1999). However, XBPs with a Be 
star companion are usually in the low-luminosity state of 
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<10 36 erg s- 1 (Bildsten et al. 1997; Stella et al. 1986), 
or a flux of < 2 x 10 -12 erg s _1 cm~ 2 at a distance of the 
SMC (60 kpc is assumed in this paper; van den Bergh 
2000). It is difficult to detect coherent pulsations from 
such a low-flux source for a nominal ASCA exposure time 
of ~ 40 ks. Therefore, we made an ASCA observation 
with an exceptionally long exposure time of ~ 177 ks. 
Within the observed region is located AX J0051.6— 7311 
= RX J0051.9— 7311, which is regarded as an XBP can- 
didate (Yokogawa et al. 2000a) , and has also been identi- 
fied with a Be star companion (Cowley et al. 1997), and 
is this the most promising target for a pulsation search. 
In this letter, we report on the discovery of coherent 
pulsations from AX J0051.6-7311 (Torii et al. 2000) with 
a very long exposure observation. We also report a long- 
term flux history based on investigating the archival data 
of ROSAT and Einstein. 

2. ASCA Observations and Data Reduction 

Three ASCA observations (hereafter obs. A1-A3) have 
covered the position of AX J0051.6— 7311. The time 
spans of these observations were 50765.179-50766.285 
(Al), 51309.604-51310.682 (A2), and 51645.986- 
51651.490 (A3) in units of Modified Julian Day (MJD). 

ASCA carries four XRTs (X-ray Telescopes, Serlemit- 
sos et al. 1995) with two GISs (Gas Imaging Spectrom- 
eters, Ohashi et al. 1996) and two SISs (Solid-state 
Imaging Spectrometers, Burke et al. 1994) on each fo- 
cal plane. In this paper, we do not refer to SIS, because 
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AX J0051. 6-7311 was outside of the SIS field of view in 
all observations. The GIS was operated in the normal 
PH mode in all observations, providing a time resolution 
of 0.0625 s (high bit rate) or 0.5 s (medium bit rate). We 
screened the GIS data by rejecting any data obtained in 
the South Atlantic Anomaly or low cut-off rigidity re- 
gions (< 4 GV). Data obtained when the satellite's ele- 
vation angle was lower than 5° were also rejected. Par- 
ticle events were removed by a rise-time discrimination 
method. After screening, the total available exposure 
times in obs. Al, A2, and A3 were ~ 43 ks, <~ 41 ks, and 
~ 177 ks, respectively. In the following analyses, we use 
the data from GIS-2 and GIS-3 simultaneously. 

3. Results 

3.1. Source Identification 

In obs. A3, a dim source was detected at ~ 13' from 
the center of the field of view (FOV). Its position was 
determined to be (00 h 51 m 38 s , -73°11'01") with an error 
of- 1'.5; thus, we designate it as AX J0051. 6-7311. This 
source was also detected in obs. Al and A2 at — 21' from 
the center of FOV, with a positional error of — l'.5-2'. 

In order to determine an X-ray counterpart, we inves- 
tigated several source catalogs of the Einstein and the 
ROSAT satellites. An Einstein source No. 25 in Wang 
and Wu (1992) and a ROSAT source RX J0051. 9-7311 
(Cowley et al. 1997; Habcrl et al. 2000) were found within 
the error region of ASCA; hence, these would be the same 
source. Cowley et al. (1997) determined the most ac- 
curate position with ROSAT HRI to be (00 h 51 m 51 s .4, 
— 73°10'38"), with a 5" error radius. Figure 1 shows a 
ROSAT HRI image around RX J0051. 9-7311 with the 
error regions of Einstein and ASCA. 

3.2. Timing Analyses 

In each observation, we collected source photons from 
a circle centered on AX J0051.6— 7311. The radius of the 
circle was set to be 2' in order to avoid contamination 
from nearby sources located > 4' away. After a barycen- 
tric correction on the photon arrival times, we performed 
an FFT (Fast Fourier Transformation) analysis. Only in 
the power spectrum of obs. A3, we detected an evident 
peak at a frequency of 0.0058 Hz. The FFT analysis was 
performed for various energy ranges to maximize the S/N 
ratio; finally, the maximum power of 59.1 was obtained 
from the photons in 2.6-8.6 keV (figure 2). The prob- 
ability to detect such a strong power at any frequency 
from random events was estimated to be ~ 8 x 10~ 8 , 
and hence the pulse detection is highly significant. We 
then performed an epoch folding search, assuming no pe- 
riod change during obs. A3. The barycentric period was 
determined to be P = 172.40 ± 0.03 s. 




Fig. 1.. ROSAT HRI image (sequence ID. rh600811n00) around 
RX J0051. 9-7311. The error circles of AX J0051. 6-7311 
(A1-A3) and Einstein source No. 25 (WW25) arc overlaid. 
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Fig. 2.. Power spectrum in obs. A3. Data points with power less 
than 8 arc omitted. An evident peak is detected at ~ 0.0058 Hz. 



Figure 3 shows the pulse profiles in the low- and high- 
energy bands (0.7-2.6 keV and 2.6-8.6 keV) and their 
intensity ratio. We clearly found the energy dependence 
of the pulse profiles; high energy photons are responsible 
for most of the pulsations. The pulsed fraction, defined as 
(pulsed flux) /(total flux) without background, is ~ 60% 
in 2.6-8.6 keV. To investigate a period change during 
obs. A3, we divided the observation span in half and per- 
formed epoch folding searches on both of the segments; 
however, no significant period change was detected. 

We also examined the aperiodic intensity variation 
during each observation, by binning the light curves with 
several time scales from — second to — hour. Neither a 
large flux variation nor any burst-like activity was found 
in each observation. Weak evidence for a random varia- 
tion of a factor — 2 was found in obs. Al and A3 with a 
time-scale of 5000-10000 s. 
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3.3. Spectral Analyses 

X-ray spectra were extracted from the same circles 
used in the timing analyses, while background spectra 
were from source-free areas near to the source regions. 

Figure 4 shows the background-subtracted and phase- 
averaged spectrum in obs. A3, which has the best statis- 
tics. Since no emission line was found in the spectra 
of obs. A1-A3, we fitted each spectrum to an absorbed 
power-law model. The spectral parameters (photon in- 
dex r and column density Ah) were found to be consis- 
tent within the statistical error throughout these obser- 
vations. Therefore, in order to obtain better statistics, 
we assumed that T and Ah are constant in these observa- 
tions, and simultaneously fitted the three spectra again. 
We obtained an acceptable fit with the best-fit parame- 
ters given in table 1. 

For obs. A3, we also extracted phase-resolved spectra 
from phases 0-0.5 ("on-pulse") and 0.5-1 ("off-pulse") 
(see figure 3). We fitted the spectra to the same model 
as given in table 1. We found that T is much smaller in 
the on-pulse phase, which is consistent with the energy- 
resolved pulse shape (figure 3). We also found a hint of 
soft excess below ~ 1 kcV only in the on-pulse spectrum. 
If we added a blackbody component to the model, the 
temperature and unabsorbed luminosity were determined 
to be 0.13 kcV and 2 x 10 35 erg s" 1 cm" 2 (0.7-10.0 keV), 
respectively, with smaller x 2 of 35.03 for 27 d.o.f.. 

4. Discussion 

The position of AX J0051.6— 7311 has previously been 
covered by 14 observations of Einstein and ROSAT. We 
used the archival data of these observations to investi- 
gate a long-term flux variation of AX J0051.6— 7311. Af- 
ter background subtraction, we derived the count rate 
of AX J0051.6— 7311 in each observation. We then con- 
verted the count rate to the flux with the PIMMS software, 
assuming that T and Ah are unchanged from the best-fit 
values for the phase-averaged spectra (table 1). The 
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Fig. 3.. Pulse profiles in the 0.7—2.6 keV (upper panel) and 
2.6—8.6 keV (middle panel) bands, and the intensity ratio be- 
tween them (lower panel). Background levels are indicated by 
broken lines. 
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Fig. 4.. Background-subtracted and phase- averaged spectrum 
from GIS2+3 in obs. A3. Crosses and a solid line indicate data 
points and the best-fit model, respectively. The best fit model 
was obtained by the spectral fitting using all the data from obs. 
Al, A2, and A3, simultaneously. 



Table 1. Results of the spectral fitting. 



Obs. ID. T JV H Flux Luminosity xV d -°- f - 

(cm -2 ) (erg s _1 cm~ 2 ) (erg s _1 ) 

— Phase-averaged — 

Al * * 6.8 x 1(T 13 2.9 x 10 35 * 

A2 * * 1.5 x 1(T 12 6.4 xlO 35 * 

A3 0.9(0.8-1.0) (< 8 x 10 20 ) 1.3 x 10~ 12 5.6 x 10 35 61.5/69 

— Phase-resolved — 

A3: on-pulse 0.6 (0.4-0.8) (< 1 x 10 21 ) 1.7 x 10~ 12 7.3 x 10 35 43.3/29 

A3: off-pulse 1.4(1.1-1.7) (< 3 x 10 21 ) 8.5 x 10~ 13 3.7 x 10 35 28.9/24 

Note — Parentheses indicate 90% confidence limits. Flux and luminosity are derived from the energy band of 0.7-10.0 keV. 
* r and Nh are linked between the observations. 
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Table 2. Flux variability of AX J0051. 6-7311. 



Obs.ID. 



Date* 
(MJD) 



Instrument 



Flux 1 " 
(erg s _1 cm" 



El. 
E2. 
E3. 
Rl. 
R2. 
R3. 
R4. 
R5. 
R6. 
R7. 
R8. 
R9. 
RIO 
Rll 
Al. 
A2. 
A3. 



43995.264 
44189.468 
44345.640 
48550.301 
48732.665 
48963.036 
49093.121 
49118.070 
49298.538 
49321.271 
49652.016 
49864.713 
50055.454 
50412.933 
50765.732 
51310.143 
51648.738 



Einstein/IPC 

Einstein/IPC 

Einstein/IPC 

ROSAT/PSPC 

ROSAT/PSPC 

ROSAT/PSPC 

ROSAT/PSPC 

ROSAT/PSPC 

ROSAT/PSPC 

ROSAT/PSPC 

ROSAT/HRI 

ROSAT/HRI 

ROSAT/HRI 

ROSAT/HRI 

ASCA/GIS 

ASCA/GIS 

ASCA/GIS 



< 1.4 x lO -1 ' 

3.4 x 10~ 13 

1.3 x 10~ 13 

8.5 x 10 -13 
7.7 x 10~ 13 

< 3.8 x 10 _1< 

< 8.6 x 10" 14 

3.4 x 10 -13 

3.4 x 10~ 13 

6.7 x 10~ 13 
8.9 x 10~ 13 

1.2 x 10~ 12 

1.8 x 10~ 12 

6.5 x 10~ 13 
6.8 x 10~ 13 
1.5 x 10~ 12 

1.3 x 10~ 12 



* Middle of the observations. 
t In the 0.7-10.0 keV band. 

results are given in table 2. We found that 

AX J0051. 6-7311 may have been brightest in the ASCA 
observations. The largest flux difference having a factor 
>50 is found between obs. R3 and obs. R10, although 
there may be some systematic error between different 
instruments (PSPC and HRI). However, even if we re- 
strict the flux comparison within the same instruments, 
we find a flux difference having a factor > 20 between 
obs. Rl and obs. R3 (PSPC). Such a large flux varia- 
tion, together with the long period (172 s), hard spec- 
trum (r = 0.9), and a Be star counterpart, suggests that 
AX J0051.6-7311 should be an XBP with a Be star com- 
panion. 

A small bump is found in the pulse profile of the lower 
energy band (figure 3; phase 0.1-0.4). The soft excess 
detected in the on-pulse phase implies existence of a pul- 
sating soft component found in some XBPs (e.g., Her 
X-l, Endo ct al. 2000; LMC X-4, Woo et al. 1996; XTE 
J0111. 2-7317, Yokogawa et al. 2000b), although better 
statistics are necessary to lead to a rigid conclusion; a 
simple power-law model is also accepted for the on-pulse 
spectrum within 95% confidence level (see table 1). 

Coherent pulsations were detected only from the data 
of obs. A3. We suggest that this is simply due to the 
much improved statistics in obs. A3, with the unusually 
long exposure time (~ 177 ks). Indeed, when we divided 
obs. A3 into four segments, each having an exposure time 
of ~ 40 ks, we could not detect significant coherent pul- 
sations from any of these segments neither by the FFT 
analysis nor by the epoch folding search. We also could 
not detect coherent pulsations from obs. R2, which has 
the best statistics among the ROSAT/Einstein observa- 



tions. Coherent pulsations have been detected from all 
of the hard X-ray sources in the SMC with a luminos- 
ity of > 10 36 erg s _1 cm~ 2 (Yokogawa et al. 2000a), 
within a typical ASCA exposure of ~ 40 ks. Pulsa- 
tions from fainter hard X-ray sources, AX J0051.6— 7311 
and 1SAX J0103. 2-7209, with luminosities of - 5 x 10 35 
erg s _1 cm~ 2 , were also detected with much longer expo- 
sures (this work and Yokogawa et al. 2000a) . These facts 
strongly suggest that a significant fraction of the uniden- 
tified hard sources in Yokogawa et al. (2000a) with a lu- 
minosity of < 10 36 erg s _1 cm~ 2 arc XBPs. Therefore, a 
long observation like obs. A3 is highly required to extend 
the pulsar population study to the low luminosity side. 

The Einstein and ROSAT data were obtained through 
the High Energy Astrophysics Science Archive Research 
Center Online Service, provided by the NASA/Goddard 
Space Flight Center. J Y. is supported by JSPS Research 
Fellowship for Young Scientists. 
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